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1.  Introduction 


In  this  p^r,  a  thenno-mechanical  problem  involving  a  finite  bar  is  considered. 
It  is  assumed  that  the  bar  initially  consists  of  two  phases,  one  of  which  was  obtained 
from  the  other  by  a  martensitic  phase  transformation.  It  is  also  assumed  that  both 
phases  of  the  bar  have  thermoelastic  constitutive  behavior.  The  two-phase  bar  is 
subject  to  dififerent  temperatures  at  each  boundary  and  to  a  boundary  traction.  A  quasi¬ 
static  thermo-mechanical  process  is  considered  and  the  temperature,  static  equilibrium 
positions,  and  quasi-static  motions  of  the  phase  boundary  are  studied  as  the  material 
coefficients,  temperature  boundary  conditions,  and  end  loading  are  varied. 

Martensitic  phase  transformations  are  diffiisionless  solid-solid  phase  transforma¬ 
tions  which  have  continuous  displacements,  with  possible  discontinuous  strains,  at 
their  phase  boundary  [2],  [S].  They  are  also  characterized  by  the  parent  and  product 
phases  having  different  crystal  structures,  or  at  least  different  crystal  oiientatioos,  and 
by  the  product  phase  having  a  shape  deformation  relative  to  the  undeformed  parent 
phase,  which  corresponds  to  an  unstressed  undeformed  state  of  that  product  phase. 
The  high-temperature  phase  of  a  material  that  can  undergo  such  a  transformation  is 
usually  referred  to  as  the  austenite  phase  of  the  material,  and  the  low-temperature 
phase  that  can  be  created  from  the  austenite  by  a  martensitic  phase  transformation 
is  usually  referred  to  as  the  martensite  phase  of  the  materiaL  Additionally,  there 
are  usually  several  variants  of  the  martensite  phase.  These  types  of  transformations 
can  be  induced  thermally  or,  in  certain  temperature  intervals,  by  the  application  of 
certain  types  of  mechanical  loadings.  In  the  absence  of  any  external  medianical 
loading,  the  temperattire  at  which  the  martensitic  phase  transformation  first  occurs  as 
the  temperature  is  lowered  is  called  the  martensitic  start  temperature  (Af,),  and  the 
temperature  at  which  the  “reverse”  martensitic  transformation  to  the  austenite  phase 


2 


first  occurs  as  the  temperature  is  raised  is  called  the  austenite  start  temperature  (A,). 
These  temperatures  are  material-dependent,  and  in  most  materials,  M,  < 


In  the  following  sections,  the  kinematics  and  the  continuum  model  that  was 
developed  in  [3]  and  is  used  here  are  discussed.  The  field  equations  for  the 
thermo-mechanical  quasi-static  process  that  is  considered  are  then  presented.  The 
constitutive  assumptions  for  each  phase  corresponding  to  the  continuum  model  that 
is  used  are  discussed.  These  constitutive  equations  are  then  linearized  about  the 
unstressed  undeformed  configurations  of  their  respective  phases.  Next,  a  kinetic 
relation  with  an  interval  of  driving  traction  conesponding  to  a  zero  phase  boundary 
velocity  is  considered.  The  linearized  forms  of  this  kinetic  relation  about  three 
different  temperatures  are  discussed.  One  of  these  temperatures  can  be  considered 
to  correspond  to  the  martensitic  start  temperature  of  the  material,  and  one  of  the 
temperatures  can  be  considered  to  correspond  to  the  austenite  start  temperature  of  the 
material.  The  boundary  value  problem  is  then  considered.  Because  the  goal  is  to 
study  the  temperature  and  quasi-static  motions  of  the  interface  and,  for  the  problem 
under  consideration,  the  displacements  are  not  needed  for  this,  only  the  temperature 
distribution  will  be  completely  solved  for.  The  temperature  at  the  interface  is  dien 
studied  as  the  parameters  of  the  problem  are  varied.  Quasi-static  motions  in  the 
neighborhoods  of  the  three  temperatures  that  were  considered  for  the  linearized  kinetic 
relation  are  then  considered.  These  quasi-static  motions  are  studied  as  the  parameters 
of  the  problem  are  varied. 
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2.  The  Kinematics 

A  one-dimensional  finite  bar  that  initially  consists  of  two  phases  is  considered.  It 
is  assumed  that  the  process  under  consideration  occurs  in  a  time  interval  F  =  [to,  ti]. 
Additionally,  for  the  problem  that  is  considered,  a  continuum  model  that  was 
developed  in  [3]  (see  also  [4])  is  used. 

Consider  a  stationaiy  reference  configuration  R  for  the  bar.  Let  x  denote  a  point 
in  R  and  let  X  be  the  length  of  the  bar  with  respect  to  this  reference  configuration. 
Considering  this,  R  can  be  expressed  as  R  =  {x /x  €  [0, X)}.  Let  y(x,t)  be  the 
suitably  smooth  and  invertible  mapping  which  maps  R  into  the  deformed  configuration 
of  the  bar  at  each  t  €  F,  with  y(x,t)  =  x  +  u(x,i)  V  (i,t)  6  R  x  F.  The  quantity 
u(x,i)  represents  the  displacement  of  a  particle  of  material  at  y  =  y(x,t)  finom  the 
point  z  6  R  at  time  t  €  F.  In  the  following,  the  two  phases  of  the  bar  will  be 
referred  to  as  phase  1  and  phase  2.  Let  x  ==  5(t)  be  the  reference  position  of  die 
phase  boundary  separating  phase  1  from  phase  2  at  time  t  €  F.  It  is  assumed  that 
particles  of  material  with  reference  points  in  R"  =  {x  /  z  £  [0,  s(t)]}  at  time  t  £  F 
are  in  phase  1,  and  it  is  assumed  that  particles  of  material  with  reference  points 
in  R'''  =  {x/z  £  [s(t),X]}  at  time  t  £  F  are  in  phase  2.  It  is  assumed  that  R" 
coincides  with  an  unstressed  undeformed  configuration  of  phase  1.  We  next  assume 
that  there  exists  a  sh^  deformation  of  phase  2  with  respect  to  R'*'  that  corresponds 
to  an  unstressed  undeformed  configuration  of  that  phase.  Let  Rf  be  the  reference 
configuration  coinciding  with  this  shape  deformation  for  all  t  £  F.  Let  zi  denote  a 
point  in  R^,  and  let  zi(z,f)  be  the  suitably  smooth  and  invertible  mapping  which 
mrqis  R"*"  into  Rj’  at  each  t  £  F,  with  Z|  =  xi(z,t)  V  z  £  R"*"  at  each  t  £  F.  For  the 
following  problem,  it  is  assumed  that  xi  is  given  by 
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x,(x,t)^x  + jmow,  (2-1) 

•(*) 

where  d(a;)  is  the  temperature  field  of  the  material  with  respect  to  the  coordinates  of  R. 
The  displacement  gradient  (transformation  strain)  corresponding  to  this  choice  of  Xi 
is  ^{6{x)),  the  Jacobian  is  J  =  1  +  7.  and  Rf  =  {xi/  xi  €  [fi(s(t),f),Zi(i,f)]}  = 
{zi/  ii  €  [s(<),£i(<)]},  where  ii(<)  =  xi(L,t).  It  is  also  required  that  7  >  -1 
so  that  reflections  are  excluded  from  (2.1).  The  mapping  which  maps  R^  into  the 
deformed  configuration  of  phase  2  at  each  t  €  P  is  represented  by  yi{xi,t),  with 
yi(xi,f)  s=  xi  +  «i(xi,f)  V  xi  e  Rj"  at  each  t  6  P. 

3.  The  Continuum  Model 

In  the  continuum  model  that  is  used  here,  the  constitutive  equations  for  each 
phase  are  defined  with  respect  to  different  reference  configurations.  More  specifically, 
the  constitutive  equations  for  phase  1  are  defined  with  respect  to  R~,  and  the 
constitutive  equations  for  phase  2  are  defined  with  respect  to  R^.  Additionally,  the 
field  equations  for  phase  1  are  expressed  with  respect  to  R~,  and  the  field  equations 
for  phase  2  are  expressed  with  respect  to  R^  (see  [3],  [4]).  The  main  advantage 
of  using  this  continuum  model  for  the  problem  under  consideration  is  that  the  field 
equations  are  in  forms  that  permit  direct  linearization.  This  is  the  case  since  the 
displacements  for  each  phase  are  measured  firom  a  reference  configuration  coinciding 
with  an  unstressed  undeformed  configuration  of  that  phase,  and  consequently,  for 
the  appropriate  boundary  and  initial  conditions,  the  displacement  gradients  can  be 
considered  infinitesimal. 
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4.  The  Field  Equations  and  Jump  Conditions 

It  is  assumed  that  the  process  under  consideration  is  thermo-mechanical  with 
no  body  forces  and  no  heat  sources  present.  It  is  further  assumed  that  the  process 
occurs  in  a  temperature  interval  Cl  =  [^g,  containing  the  temperatures  M,  and  A,. 
In  the  following,  a  quasi-static  process  is  considered.  In  particular,  a  process  is 
considered  where  time  is  treated  as  just  a  parameter  in  all  quantities  except  for  s(t).^ 
The  dependence  of  s(t)  on  time  is  such  that  (n  >  1)  and  S(f)  =  are 

assumed  to  be  negligible  for  all  t  €  F,  and  are  therefore  taken  to  be  equal  to  zero. 
For  this  type  of  process,  the  only  time  derivative  terms  in  the  boundary  value  problem 
that  are  not  identically  equal  to  zero  are  those  that  are  linear  in  i(t)  =  Such  a 
process  corresponds  to  a  process  where  |i(f)|  «  1  and  the  acceleration  of  the  phase 
boundary  is  negligible  for  all  f  €  F. 

The  general  held  equations  using  the  type  of  continuum  model  described  in 
previous  section  were  derived  and  discussed  in  [3].  For  the  type  of  process  under 
consideration,  the  balance  of  linear  momentum,  the  balance  of  energy  (the  first  law 
of  thermodynamics),  and  the  entropy  inequality  (the  second  law  of  thermodynamics) 
for  phase  1  are 


(4.1) 


_ 

^  Becsoie  of  this,  in  the  following,  the  dependence  of  all  quantities,  except  for  s(()  and  xi(z,  (),  on  time  will 
not  be  displayed  in  the  notation.  Also,  strictly  freaking,  such  a  process  is  not  a  true  quasi-static  process  since 
time  is  not  just  a  parameter  in  all  variables  of  the  problem,  and  thmfore  the  set  of  all  solutions  as  time  is  varied 
doea  not  consist  of  only  static  equilibrium  solutions. 
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respectively,  V  x  €  (0,  ^(t))  at  each  t  6  F,  where  <t(x)  is  the  nominal  stress  with 
respect  to  R~  for  phase  1,  g(x)  is  the  nominal  heat  flux  with  respect  to  R~  for  phase 
1,  and,  as  mentioned  previously,  6(x)  is  tte  temperature  field  of  the  material  with 
respect  to  the  coordinates  of  R. 

Similarly,  for  the  type  of  process  under  consideration,  the  balance  of  linear 
momentum,  the  balance  of  energy,  and  the  entropy  inequality  for  phase  2  are 


dqt 

dxt 


=  0, 


(4.2) 


jI 

dxi 


<0, 


respectively,  Vzi  €  £i(f))  at  each  f  €  F,  where  ert(xi)  is  the  stress  with  respect 

to  Rf  for  phase  2,  9i(zi)  is  the  heat  flux  with  respect  to  R|  for  phase  2,  and  0i(xt) 
is  the  temperature  field  of  phase  2  with  respect  to  the  coordinates  of  R^. 

In  addition  to  continuity  of  displacements  at  the  phase  boundary,  we  also  require 
that  the  temperature  be  continuous.  At  x  =  s(i),  these  conditions  can  be  expressed 
as 


0+(s(<))  =  e-(s(t)), 


(4.3) 


respectively,  at  each  t  €  F,’  where  the  and  superscripts  denote  the  limiting 

*  Note  that  the  original  fonn  of  the  continuity  of  (Usplacement  oooditkn  yf  =  y~  redacea  to  (4.3)i  for  the  ii 
given  by  (2.1). 
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values  of  their  respective  quantities  from  the  “4-"  and  sides  of  i  =  ^(t), 
respectively,  at  time  t  6  F. 

As  mentioned  in  the  Introduction,  the  strains  at  the  phase  boundary  separating 
two  phases  involved  in  a  martensitic  phase  transformation  may  be  discontinuous. 
Considering  this,  we  require  that  the  first  and  second  derivatives  of  y  be  piecewise 
continuous  on  R,  with  possible  discontinuities  occurring  only  at  z  =  s{t).  In  adcUtion 
to  this,  we  require  that  the  stress,  entropy,  and  heat  flux  be  piecewise  continuous,  with 
possible  discontinuities  occurring  only  at  z  =  s(i).  The  jump  conditions  at  z  ==  s(t) 
that  are  equivalent  to  the  balance  of  linear  momentum,  the  balance  of  eneigy,  and  the 
entropy  inequality  are  given  by 

-<T~  -  0, 


-  V")  +  f}s  +  {qt  -  q~)  =  0,  (4.4) 

fa  >  0, 

respectively,  at  each  t  €  F,  where  ?;(z)  is  the  entropy  of  phase  1  with  respect  to  R~, 
r7i(zi)  is  the  entropy  of  phase  2  with  respect  to  R^,  and 

/  =  J0,-*--^--<7~(7  +  7+7  +  7+-7~)  (4J) 

is  the  driving  traction  acting  at  the  interface  for  the  quasi-static  process  under 
consideration.^  We  also  have  at  the  phase  boundary  a  kinetic  relation  [1]  relating 
i  to  at  least  /.  This  kinetic  relation  is  a  constitutive  relation  and  will  be  discussed 
in  Section  6. 

*  See  [1]  for  the  fomu  of  the  field  equatioof  and  jwsp  conditiona  that  were  pieaenled  in  dua  aection  that 
coneapond  to  a  contimmm  model  that  uaea  a  aingle  elaitic  potential  and  only  one  lefieienoe  ooofignntion. 
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We  must  also  specify  the  boundaiy  conditions  at  x  =  0  and  xi  =  Li(t).  The 
boundary  conditions  for  Equations  (4.1)]  and  (4.2)i  are 

u(0)  =  0,  ai(Li(t))  =  (To  5  (4.6) 

respectively,  and  the  boundary  conditions  for  the  temperatures  6(x)  and  0i(x,)  are 

<?(0)  =  d-,  9y(Li(t))  =  9: ,  (4.7) 


respectively. 


5.  The  Constitutive  Assumptions 

It  is  assumed  that  both  phase  1  and  phase  2  are  homogeneous  thermoelastic 
phases.  In  particular,  for  phase  1  it  is  assumed  that  there  exists  a  Helmholtz  free 
energy  potential 


^  =  ^(7,0), 


such  that 


7  = 


dxj; 

W 


(5.1) 


(5.2) 


where  y  =  and  it  is  assumed  that  the  heat  flux  response  function  for  phase  1 
has  the  form 

s)- 

Similarly,  for  phase  2  it  is  assumed  that  there  exists  a  Helmholtz  free  energy  potential 


(5.4) 
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such  that 


dyi>i 


drift 

dOy' 


(55) 


where  71  =  and  it  is  assumed  that  the  heat  flux  response  function  for  phase  2 
has  the  form^ 


We  note  that  all  of  the  quantities  for  phase  1  given  by  (5.1H3.3)  are  defined  with 
respect  to  the  reference  configuration  R~,  and  all  of  the  quantities  for  phase  2  given 
by  (S.4H5.6)  are  defined  with  respect  to  the  reference  configuration  Rf . 

We  next  assume  that  the  thermo-mechanical  process  under  consideration  talces 
place  in  a  neighborhood  of  a  temperature  60  €  il.  More  precisely,  letting  0  ^ 
(0  —  0q)  and  01  =  ((>1  —  5o)»  we  assume  that  1 0 1  «  1  Vx  €  [0, s(/)]  and  l^il  «  1 
Vzi  €  [s(<),  Ii(t)],  at  each  t  €  F.  Additionally,  we  assume  that  for  a  given  ^o,  the 
dependence  of  7  on  d  is  at  least  second-order;  i.e.  we  assume 

7  =  70  +  0(^’),  (5.7) 

where  70  =  7(^0)  >  —1-  For  this  assumption,  the  first-order  approximation  of  zi  is 

xi(x,  t)  =  X  -f-  7o(x  -  s(t)),  (5,8) 

and  the  first-order  {q)proximation  of  L\{t)  is  L\{t)  =  L  y(t{L  —  s(i)).  We 
assume  that  phase  1  is  unstressed  at  (7,0)  =  (O,0o).  and  phase  2  is  unstressed 


’  The  constitutive  equations  conesponding  to  a  thennoelastic  material  defined  with  lespect  to  a  single  fixed 
leCerenoe  confignntioa  are  well  known  and  can  be  found  in,  e.g„  [1]. 


at  (7ii^i)  =  (0,^o)-  We  also  assume  that  j^l  «  1  Va:  e  [0,j(0)  and  |7i|  «  1 
Vzi  €  {s{t),  >Lj(0]>  at  each  t  €  F.  Considering  this,  i?  for  phase  1  can  be  written  as 

^  =  +  0(d^,7d,7^),  (5.9) 

where  V'*  =  0(0,  ^o)  and  i/’  =  t](Q,0q)  =  -|f|(o,sp)-  Similarly,  0i  for  phase  2 
can  be  written  as 

01  =  0r  -  n‘Bi  +  o(^?,7i^„7j*^,  (5.10) 

where  0*  =  0i(O,do)  and  rf’  =  >?i(0,6o)  =  —  ^  •  We  next  assume  that  the 

*  l(0,flo) 

heat  flux  q  for  phase  1  is  zero  when  ||  =  0,  and  the  heat  flux  91  for  phase  2  is 
zero  when  sj-  =  0.  We  additionally  assume  that  |s|  «  1  V  z  e  [0,s(<))  and 
«  1 V  €  (a(t),  li(t)],  at  each  t  €  F.  For  these  assumptions,  the  first-<»der 
approximation  of  q  for  phase  1  is 

«  =  (5.11) 

ax 

where  k  =  [dq/d{^)]  ^  and  the  flrst-K)rder  approximation  of  91  for  phase  2  is 

9i  =  fc.^,  (5.12) 

(ix\ 

where  kt  =  ^  *  Using  (5.9),  (5.10),  and  (4.3)2,  the  first-order 

approximation  of  the  driving  traction  given  by  (4.5)  is 

/  =  J-0:  -  0-  -  (JV  -  rj-)0^  -  a-70 ,  (5.13) 

where  a~  is  given  by  the  linearized  form  of  (5.2)i  and  J*  =  I  +  7(^0)  =  1  +  7o 
We  must  also  specify  the  kinetic  relation,  which  relates  i  and  at  least  /.  This  will 
be  discussed  in  the  following  section. 

*  Note  that  Equatkns  (5.1 1)  and  (5.12)  both  have  the  foim  of  Fouiier’a  Law. 

^  Note  that  became  both  phases  are  homogeneom,  =  i>*,  n~(0,So)  =  e*.  = 

tfr,*,  and  t»,^(0,do)  =  n,*. 
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6.  The  Kinetic  Relation 

It  is  assumed  that  the  thermoelastic  material  under  consideration  has  a  kinetic 
relation  of  the  form 


(6.1) 

where  0  =  0(3).  Because  of  the  entropy  inequality  (4.4)3,  this  kinetic  relation  must 
be  such  that 


mf,0)>O. 


(6.2) 


We  additionally  assume  that  at  each  temperature  0(s)  €  H  the  kinetic  relation  given 
by  (6.1)  has  the  form 


*(/.«)  1 


V 


>0, 
=  0, 
<0, 


/  >  fm 

MB)  <f<  MB)  , 
/  <  MB) 


(6.3) 


where  /i(d(s))  and  MB(^))  depend  on  the  material  and  are  such  that  /i  <  0  and  > 
0  V  d(s)  €  fl.  We  note  that  this  kinetic  relation  satisfies  (6.2).  A  kinetic  relation 
that  has  the  form  given  by  (6.3)  and  is  linear  in  /  is 


/  >  MB) 

MB)  <f<  MB)  ,  (6.4) 

/  <  MB) 

where  J^i(^(s))  >  0  and  »'j(^(j))  >  0  V  9{3)  €  ft.  We  note  that  because 
(6.4)  is  linear  in  /,  it  is  probably  most  appropriate  only  for  values  of  /  such 
that  /i  <  /  <  /2,  !(/  -  /a)//!  «  1,  or  |(/  -  fi)/f\  «  1.  For  the 


f  1 


MB) 


{ 


I  MB) 


{/  -  MB)h 

0, 

{/  -  MB)}, 
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foUowing,  let  B(do)  =  {(7■>7^^)/  |7~I  «  1,  l7ri  «  1, 1^1  «  l}  and  A(0o)  = 
{(/(7~»7r>^)>^)/  (7"»7i*‘>^)  €  B(6o)}.  for  a  given  €  fi.*  Three  cases  will  be 
considered. 


Case  r  Consider  first  the  case  where  the  temperature  and  the  material 

are  such  that  /(0,0,^o)  =  fi(^o)  =  M^o)  =  0’  Additionally,  assume  that 
€  B(0o)  at  each  t  e  P.  For  this  case,  =0  and  ^  =  0, 

since  fi(0o)  =  M^o)  =  0  and  h  is  required  that  /i  <  0  and  /s  >  0  V  0(s)  € 
Thus,  for  this  case,  /i  =  O(d^),  f2  =  0(0^),  and  the  first-order  approximation  of 
(6.4)  is 


[  -  (J-T,-  -  +  a-7o},  -( -  r)-)0  ~  a-70  >  0 

^  ~  1  ’ 

[  +  -  f'lc  <  0 

(6^) 

(7‘'i7i'‘i^)  €  B(0o).  where  i/j  =  ui{0o)  >  0  and  i/j  =  >'»(0o)  >  0- 

Case  II:  (Consider  the  case  where  the  temperature  ^  and  the  material  are  such 

that  /(0,0,^o)  =  /i(^o)  and  /i(^o)  **  not  small,  in  the  sense  that  \f  -  fi\  <  fi 

V  (/,d)  €  A(do)'  If  phase  1  represents  the  austenite  and  phase  2  represents  the 

martensite,  the  temperature  0o  for  this  case  might  be  considered  to  represent  the 

martensitic  start  temperature  of  the  material.  We  also  assume  that  (7*,  'y^,0)  e  B(do) 

at  each  t  €  T.  For  this  case,  the  portion  of  (6.4)  that  corresponds  to  (/,  0)  €  A(do)  is 

r  0,  /  >  m 


HM={  I 


/</iW 


Substituting  (5.13)  and 


/i  =  /i(W  +  /<i«  +  0(*'), 


'  Notethatalthough/depeii(UaaS(«).itcaiibevanedindepeiideDUyofS(i)unoeUilaodepeiidaoo7  andyi*^. 
*  Note  that  /(O.O.So)  = 
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where  ^  ,  into  (6.6)  and  retaining  only  first-order  terms  yields 


0, 


$  = 


(6.8) 


V(7-,7^^)GB(^o). 


Case  ni:  We  next  consider  the  case  where  the  temperature  C  and  the  material  are 
such  that  /(O,O,0o)  =  M^o)  and  f^ido)  is  not  small,  in  the  sense  that  \f  -  /2I  <  /? 
V  if, 6)  e  A(do)<  If  phase  1  represents  the  austenite  and  phase  2  represents  the 
martensite,  the  temperature  6q  for  this  case  might  be  considered  to  represent  the 
austenite  start  temperature  of  the  material.  We  also  assume  that  (7~,7i''',  B)  €  B(^) 
at  each  ^  €  F.  For  this  case,  the  portion  of  (6.4)  that  corresponds  to  (/,  0)  6  A(^o)  is 


^im  = 


if  -  f7(e)), 

0, 


/  >  M0) 
f  <  f2i0) 


(6.9) 


Substituting  (5.13)  and 


/r  =  /2(^o)  +  f^70  +  0(0* ), 


(6.10) 


where  /i,  =  ^  ,  into  (6.9)  and  retaining  only  first-order  terms  yields 

^<1 


0, 


-(j"fli-fl'‘)0-<r~7o  <  f*70 

(6.11) 


V  (7-,7^<?)  €  B(0o).»° 

Note  that  cbeie  m  other  cuts  conceniijis  the  iioeariaed  form  of  (6.4)  that  can  be  oonsideied.  For  eaampte, 
one  can  consider  the  case  wheie  and/or  are  small  and  nonzero. 
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Another  method  for  obtaining  a  linear  approximation  of  the  nonzero  portion  of 
a  kinetic  relation  of  the  form  (6.1)  is  to  simply  expand  it  in  a  Taylor  series  about 
{f,6)  =  {fo,do),  where  /o  =  /(0,0,^o)-  For  the  moment,  assume  that  (6.1)  and  its 
first  derivatives  are  continuous  at  (/,  6)  =  (/o,  ^o)-  For  the  case  where  #(/o,  ^o)  =  0. 
the  Taylor  series  expansion  of  (6.1)  about  (f,6)  =  {fo,6o)  is 

$  =  «.(/-  /o)  +  M  +  O  ((/  -  /o)^  (/  -  fo)e,  ,  (6.12) 

where  a,  =  |jl(/o,flo)  and  Oj  =  |||(/o,tfo)-”  Substituting  (5.13)  into  (6.12)  and 
retaining  only  first-order  terms  yields 

$=  (6.13) 

For  the  case  where  aj  is  not  continuous  at  /o,  a]  is  taken  to  be  the  limiting  value 
of  ly  as  /  approaches  /o  from  the  side  of  /  corresponding  to  a  nonzero  portion  of 
$.  Comparing  (63)  with  (6.13),  we  can  conclude  that  Q|  —  ^  and  a,  —  0,  where 
1/  s  i/t  or  t/  =  1/3,  depending  on  which  branch  of  (6.5)  is  considered.  Comparing 
(6.8)  and  (6.11)  with  (6.13),  we  can  conclude  that  ai  =  ^  and  aj  =  -a,  where 
u  —  Ui,  fi  =  Hi  or  u  =  Ui,  /i  =  /i],  depending  on  the  case  considered. 

7.  The  Solution  of  the  Boundary  Value  Problem 

From  Equation  (4.2)|  and  the  boundary  condition  given  by  (4.6)2,  we  can  conclude 
that 

<7i(xi)  =  <7o  (7.1) 

for  s{t)  <  xi  <  Li{t),  at  each  t  €  F.  From  Equation  (4.1)i,  the  jump  condition 
(4.4)i,  and  (7.1)  given  above,  we  obtain 

a(x)  =  <70 ,  (7.2) 

"  Note  that,  although  /  -  /o  it  not  small  for  a  relatively  large  range  of  values  of  (y~,  0)  €  f  -  fa 

does  go  to  aero  at  the  tame  rate  u  y~,  7^,  aQdS(a)  go  to  reia 


for  0  <  X  <  s{t),  at  each  t  eT.  Since  the  objective  is  to  calculate  the  temperature 
at  the  interface  and  the  corresponding  kinetic  relation,  the  displacements  will  not  be 
calculated. 


Using  (5.11),  Equation  (4.1)2  for  phase  1  becomes 


(7.3) 


Vi  6  (0,a(t)),  at  each  t  eT.  Using  (5.11)  and  (7.3)  given  above.  Equation  (4.1)3 
becomes 

/ 1 

(7.4) 


-(«) 


Vi  €  (0,s(t)),  at  each  f  e  F.  Additionally,  the  boundary  condition  (4.7)i  in  terms 
of  6  is 

^‘(0)  =  6,  (7.5) 

where  6  =  9*- 9q.  From  (7.4),  we  obtain  the  well  known  result  that  1:  >  0.  Assuming 
that  ib  >  0,  the  solution  of  (7.3)  that  satisfies  the  boundary  condition  (7.5)  is 

^(i)  =  Cn  +  9,  (7.6) 

for  0  <  I  <  s(t),  at  each  t  €  F. 

Using  (5.12),  Equation  (4.2)2  for  phase  2  becomes 


=  0, 


(7.7) 


Vii  e  (a(t),Li(<)),  at  each  t  €  F.  Using  (5.12)  and  (7.7)  given  above.  Equation 


(4.2)3  yields  ki  >  0.  The  boundary  condition  (4.7)2  in  terms  of  (9i  is 


l(Li(t))  =  ^1,  (7.8) 

where  -  $q.  Assuming  that  ib|  >  0,  the  solution  of  (7.7)  that  satisfies  the 

boundary  condition  (7.8)  is 

^1(2^1)  ~  ^3(^1  —  Li(t))  +  ,  (7.9) 

for  s(t)  <  xi  <  Li(t),  at  each  <  e  F.  From  the  continuity  of  temperature  jump 
condition,  we  can  solve  for  Ci  in  (7.6)  in  terms  of  C3  in  (7.9)  and  obtain 

^(x)  =  {-C3(l  +  7o)(L (7.10) 

s 

for  0  <  X  <  s(<),  at  each  t  6  F. 

For  the  calculation  of  the  first-order  ^proximation  of  the  energy  jump  condition, 
a  kinetic  relation  which  can  have  the  forms  of  the  kinetic  relations  given  by  (6.5), 
(6.8),  and  (6.11)  will  be  used.  In  particular,  we  consider  the  kinetic  relation  which 
has  the  form 

i  =  -/•(*))  •  ail) 

We  next  assume  that  the  temperature  and  the  material  are  such  that  /(0, 0,  ^o)  = 
/*(^o)*  Considering  this  and  substituting  (5.13)  and 

/■  =  m)  +  »*9+o(«»),  am 

where  ft  =  ^  ,  into  (7.1 1)  and  retaining  only  fiiat-order  terms  yields 
00 

i  =  [m  +  (J'nr  -  >)•)]»  +  <r'7o},  (7.13) 

€  B(0o),  where  a"  =  (7o  fof  the  problem  considered  here.  This  kinetic 
relation  will  be  used  for  the  boundary  value  problem  under  consideration.  Using 


(5.13),  (7.13),  (5.11),  and  (5.12),  the  iirst-order  approximation  of  the  energy  jump 
condition  given  by  (4.4)2  is 


-»?■■)  +/o}{[m+  (/*»?r-»7*)]^  +  <yo7o} 


(7.14) 


The  constant  C3  in  Equations  (7.9)  and  (7.10)  will  be  determined  from  this  equation. 
We  next  define  the  nondimensional  quantities 


(To  = 


II 

II 

«•  - 

J-k 

(7.15) 

<To70 

uk 

^0  [n  +  ( L  ( -  £•)  n+ 


where  is  the  internal  energy  of  phase  2  evaluated  at  (71,0,)  =  (0,^o).  e*  is 
the  internal  energy  of  phase  1  evaluated  at  (7, 9)  =  (0, 9o),  and  ( -  e*)  = 

(j'Vr  -  +  0o[j*Vi  -  ^*).  Substituting  (7.9)  and  (7.10)  into  (7.14),  using  the 

nondimensional  quantities  given  by  (7.15),  and  then  solving  for  C3  yields 


^3- 


id +  £  + 


(7.16) 


Substituting  (7.16)  into  (7.9)  or  (7.10)  jdelds  the  following  for  the  tenq)eratuie  at  the 
interface: 
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~  S)+i>  k  §*S 
(i/  +  S)(l  —  s)+i/kS 


(7.17) 


8.  The  Temperature  at  the  Interface 

The  temperature  at  the  interface  given  by  (7.17)  for  the  three  cases  described 
in  Section  6  will  be  studied  as  the  material  coefficients  i>  and  k,  the  temperature 
boundary  conditions  6*  and  d*.  the  end  load  ao,  and  the  position  i  of  the  phase 
boundary  are  varied.  The  temperature  at  the  interface  for  each  of  these  three  cases 
can  be  obtained  from  (7.17)  by  substituting  the  appropriate  values  of  u  and  fi.  For 
Cases  II  and  m,  when  /  is  in  an  interval  such  that  i  =  0,  the  temperature  at  the 
interface  can  be  obtained  from  (7.17)  by  letting  v  oo. 

Plots  of  d(a)  vs  a  for  several  values  of  v  are  presented  in  Figures  1  and  2.  In  these 
figures,  and  in  the  rest  of  the  hgures  in  this  section,  the  values  of  ^(s), 
that  are  shown  should  all  be  multiplied  by  some  small  number,  e.g.  10~^.  We  note 
that  such  a  multiplicative  factor  can  be  divided  out  in  (7.17).  The  limiting  values  of 
6(3)  as  i>  — »  0  and  x/  — ►  00  are 


lim  ^(i)  =  — ctq  , 


lim  0(8)  = 

v-*oo 


0*(1  -3)  +  9\k3 
I  —  5+ ka 


(8.1) 


respectively.  An  easy  calculation  shows  that  as  1/  is  varied  0  only  as 

v  -4  00,  and  that  =  0.  Therefore,  as  i/  is  varied,  6(a)  (at  a 

given  1)  nwnotonically  increases  or  decreases,  depending  on  the  values  of  the  other 
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parameters. Also,  note  that  (8.1)2  is  the  temperature  at  the  interface  for  Cases  n 
and  m  when  /  is  in  an  interval  such  that  i  =  0. 

Plots  of  ^(s)  vs  3  for  several  values  of  ib  are  presented  in  Figure  3.  The  limiting 
values  of  ${s)  as  ib  -+  0  and  ib  -»  00  are 


fc— 0  S  +  1/ 

(8.2) 

lim  #(s)  =  6; , 

fc-*oo 


respectively.  An  easy  calculation  shows  that  as  ib  is  varied  »  0  only  as  ib  — »  00, 
and  that  limi;_oo  =  0.  Therefore,  as  ib  is  varied,  ^(5)  monotonically  increases 
or  decreases,  depending  on  the  values  of  the  other  parameters. 

Plots  of  0(S)  vs  s  for  several  values  of  ^0  are  presented  in  Figures  4  and  5.  The 
limiting  value  of  d(s)  as  ao  — »  0  can  easily  be  observed  from  (7.17),  and  the  limiting 
value  of  0(5)  as  ag  -+  00  is  00. 

The  values  of  0(5)  as  5  — »  0  and  5  — »  1  are 


lim 0(5)  =  0*,  lim 0(5)  =  0,  , 


(8.3) 


respectively.  These  results  should  be  expected.  An  extremum  of  0(5)  can,  of  course, 
be  determined  from 


<f0(5) 

ds 


=  0. 


(8.4) 


For  a  given  set  of  material  coefficients  and  boundary  conditions,  the  values  of  5  that 
satisfy  (8.4)  are 


Whether  0{$)  monotonically  increases  or  decreases  can.  of  course,  be  determined  by  comparinf  (8.1))  and 
(8.1)2  *  8>ven  set  of  parameter  values. 
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3  = 


ao  + 


cro(l-it) 


ao  + 


vk{9*-  -  (ao  + 

^  —  k6i  j  ao^l  -  +  ^*  —  kB^ 


(83) 


Note  that  for  a  given  set  of  parameters,  a  given  by  (8.5)  is  a  physically  meaninghil 
solution  only  if  it  is  real-valued  and  is  such  that  0  <  I  <  1. 


9.  A  Quasi-Static  Process 

In  this  section,  a  quasi-static  process  will  be  considered  for  the  three  cases 
described  in  Section  6.  The  direction  of  motion  of  the  phase  boundary  will 
be  considered  and  the  static  equilibrium  positions  of  the  phase  boundary  will  be 
determined.  For  the  following  cases,  let 


i  ds(i) 


(9.1) 


Before  we  consider  the  three  cases,  however,  we  will  first  consider  the  signs  of  some 
of  the  variables  that  have  the  same  signs  for  all  of  the  cases  considered.  As  is  usually 
done,  we  assume  that  the  dependence  of  the  Gibbs  firee  energy  g*  for  an  unstressed 
phase  1  and  the  Gibbs  free  energy  for  an  unstressed  phase  2  (with  respect  to  R^) 
on  the  temperature  is  qualitatively  like  that  presented  in  Figure  6  (see  [2]).'^  Because 
the  Gibbs  free  energy  of  an  unstressed  material  is  equal  to  the  Helmholtz  firee  energy 
of  that  unstressed  material,  we  have  —  g*  and  for  all  temperatures 

The  Gibbs  free  raergy  of  the  phase  1  can  be  defined  asp  =  t^-0^Tandthe  Gibbs  free  eneigy  of  phase  2  can 
be  defined  as  pi  =  sfti  -  a\i\. 
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where  the  undeformed  configurations  of  both  phases  are  unstressed.  From  this,  we 
can  con  lade  that  >7’  >  J’t/*  V  d  €  fi  for  the  case  under  consideration.  We  also 
always  have  L  >  Q,  ^  0i  J  >  0,  i  >  0,  >  0,  and  u  >  0.  From  this,  we 

can  conclude  that  ib  >  0.  We  can  also  conclude  from  this,  and  the  kinetic  relations 
presented  in  Section  6,  that  for  a  given  <7o,  an  increase  in  temperature  will  result 
in  an  increase  in  /,  and  a  decrease  in  temperature  will  result  in  a  decrease  in  /. 
Therefore,  if  phase  1  represents  the  austenite  and  phase  2  represents  the  maitensite, 
which  is  what  corresponds  to  Figure  6,  an  increase  in  temperature  favors  the  austenite, 
and  a  decrease  in  temperature  favors  the  martensite.  This  corresponds  exactly  with 
experimental  observations  [2],  [5]. 

Case  I:  From  the  description  of  this  case  given  in  Section  6,  we  have  J'V’’  -  = 

0  and  /i  =  0.  We  next  make  the  additional  assumption  that  1/1  =  1/,  =  v,  and  we  let 

=  »'L=o »  ^0  =  •  (^-2) 

For  this  case,  the  sign  of  i  is  opposite  the  sign  of  i,  the  sign  of  ctq  is  opposite 
the  sign  of  ao7o,  i>'  >  0,  and  the  kinetic  relation  given  by  (6.S)  in  terms  of  the 
nondimensional  variables  is 

(9J) 

where  ^  is  the  temperature  at  the  interface  given  by  (7.17)  with  (P,ao)  = 

For  the  quasi-static  process  under  consideration,  static  equilibrium  occurs  when  i  =  0. 
From  (9.3),  we  can  conclude  that  this  occurs  when 

9(3)=  (9  Ay* 

In  the  remainder  of  this  lection.  the  temperature  at  the  interface  may  lometimet  be  referred  to  at  aimply  tf(a). 
with  the  undeistandins  that  it  corresponds  to  the  case  under  consideration. 
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Using  (7.17)  in  (9.4),  the  static  equilibrium  position  of  the  phase  boundary  for  a  given 
set  of  material  coefficients  and  boundary  conditions  will  occur  at 


3  = 


6*  —  k6*  + 


(9.5) 


with  the  requirement  that  0  <  3  <  1.  We  note  that  in  Figures  1  and  2,  where  plots  of 
6(1)  vs  3  for  several  values  of  t/  are  presented,  the  static  equilibrium  positions  given 
by  (9.5)  occur  where  all  of  the  temperature  curves  intersect.  This  is  because  the  static 
equilibrium  position  of  the  phase  boundary  given  by  (9.5)  is  independent  of  t/.  From 
(9.3)  and  the  fact  that  the  sign  of  3  is  opposite  the  sign  of  i,  we  can  conclude  that 
3  >  0  for  6(3)  >  -  ctq  ,  and  i  <  0  for  0(s)  <  —  Because  of  this,  the  static 
equilibrium  positions  in  Figures  1  and  4  are  unstable,  in  the  sense  that  the  phase 
boundary  will  move  away  from  these  static  equilibrium  positions  for  any  perturbation 
away  from  them,  and  the  static  equilibrium  positions  in  Figures  2,  3,  and  5  are  stable. 
In  Hgures  4  and  5,  some  values  of  are  used  which  result  in  the  static  equilibrium 
position  determined  by  (9.5)  being  outside  of  the  interval  (0, 1).  Consequently,  the 
static  equilibrium  states  of  the  bar  for  these  values  of  consist  of  only  one  phase. 

The  case  where  j*r/*  =  r/*  in  addition  to  =  0*  can  also  be  considered 
hete.'^  In  particular,  for  this  case,  the  first-order  ^proximation  of  the  energy  jump 
condition  reduces  to  the  continuity  of  heat  fluxes  across  the  phase  boundary.  From 
(7.14),  we  can  conclude  that  the  temperature  solution  for  this  case  is  the  same  as  the 
temperature  solution  for  the  case  where  u  oo.  Thus,  the  temperature  and  quasi¬ 
static  motions  for  the  case  where  =  0’  and  J*»;*  =  t;*  are  the  same  as  they  are 
for  the  cases  presented  here  for  Case  1  where  v  oo. 

These  astumpiioa  might  be  most  spproprisie  to  the  esse  when  phase  t  and  phase  2  repieseat  two  different 
valiants  of  the  same  martensite. 
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Case  U:  For  this  case,  let 


^1  = 


(9.6) 


From  the  description  of  this  case  given  in  Section  6,  we  have  <  V**.  If  we  also 
make  the  assumption  that  |/i|  <  {J’rj*  —  r/’l,  the  sign  of  3  will  be  opposite  that  of  i, 
the  sign  of  ay  will  be  opposite  that  of  ao'yo,  and  i>i  >  0.  The  kinetic  relation  given 
by  (6.8)  in  terms  of  the  nondimensional  variables  is 


(9.7) 


.  1 

(«® + «.) 

+  ai) 

where  is  the  temperature  at  the  interface  given  by  (7.17)  with  (i>,  ao)  —  (i'n  o^i)» 
and  is  given  by  (8.1)2,  which  corresponds  to  (7.17)  with  x/  00.  In  Figures 
7  and  8,  plots  of  &(3)  vs  3  for  u  =  0.1  and  — »  00  are  presented.  The  case 
where  r/  =  0.1  represents  with  i/y  =  0.1,  and  the  case  where  1/  — ►  00  represents 
Quasi-static  processes  will  be  considered  for  the  cases  presented  in  these  two 
figures.  The  position  of  the  phase  boundary  corresponding  to  0(3)  =  —  ao  is  given 
by  (9.5)  with  00  =  ^1.  This  static  equilibrium  position  also  corresponds  to  the  point 
where  0^^^  and  0^^^  intersect,  since  it  is  independent  of  i/.  When  0(3(to))  <  -^1, 
the  temperature  at  the  interface  is  given  by  and  when  0(3 (to))  >  -ay,  the 
temperature  at  the  interface  is  given  by  O^^K  From  (9.7),  we  can  conclude  that  a 
phase  boundary  with  an  initial  position  of  S(fo)  such  that  0(5(to))  >  —^1  will 
stay  at  that  position,  and  a  phase  transformation  will  not  occur.  However,  a  phase 
boundary  with  an  initial  position  such  that  0(3(io))  <  —ay  will  have  a  i  <  0,  and 
consequently  phase  1  will  be  converted  into  phase  2. 
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Case  ni:  For  this  case,  let 


From  the  description  of  this  case  given  in  Section  6,  we  have  .  Also,  for  die 

following,  we  assume  that  \n\  <  |  J*t7*  —  rf^\.  For  this  case,  the  sign  of  u  and  the  sign 
of  i(f)  in  relation  to  S  depend  on  the  relative  values  of  (Jip*  -  V’*)  and  doiJql  ~  7*)- 
Thus,  we  cannot  make  such  definite  conclusions  concerning  the  sign  of  u  and  the  sign 
of  i(t)  in  relation  to  s  for  this  case,  as  we  were  able  to  do  for  the  two  cases  discussed 
previously.  The  kinetic  relation  given  by  (6.11)  in  terms  of  the  nondimensional 
variables  is 


(9.9) 


where  is  the  temperature  at  the  interface  given  by  (8.1)2,  which  corresponds  to 
(7.17)  with  1/  -»  00,  and  is  the  temperature  at  the  interface  given  by  (7.17)  with 
(i/, ^q)  =  (i/,,  a2).  For  this  case,  we  will  consider  a  quasi-static  process  for  the  cases 
presented  in  Figures  7  and  8.  The  position  of  the  phase  boundary  corresponding 
to  B(s)  =  -  cTj  is  given  by  (9.5)  with  a'  =  ^2.  This  static  equilibrium  position 
also  corresponds  to  the  point  where  and  intersect,  since  it  is  independent  of 
P.  When  0(3 (to))  <  -ai,  the  temperature  at  the  interface  is  given  by  and 
when  d(5(fo))  >  -  ^i,  the  temperature  at  the  interface  is  given  by  From  (9.9), 
we  can  conclude  that  a  phase  boundary  with  an  initial  position  of  5(<o)  such  that 
0(s(fo))  <  -  ^2  will  stay  at  that  position,  and  a  phase  transformation  will  not  occur. 
However,  a  phase  boundary  with  an  initial  position  such  that  tf(s(fo))  >  will 
have  a  i  >  0,  and  consequently  phase  2  will  be  converted  into  phase  1. 
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